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Abstract 



By combining the inputs from the neutrinoless double beta decay and 

the fits of cosmological models of dark matter with solar and atmospheric 

neutrino data, we obtain constraints on two of the mixing angles of Majorana 

neutrinos, which become stronger when coupled with the reactor neutrino 

data. These constraints are strong enough to rule out Majorana neutrinos if 

the small angle solution of solar neutrino puzzle is borne out. 
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It is well known [[T]|| that unless neutrinos are very massive and nonrelativistic, or 
interact through both left and right-handed currents, experimental data on neutrino-induced 
reactions cannot distinguish between Dirac and Majorana neutrinos. Neutrinoless double 
beta decay remains as the only feasible tool to probe this question. Although experiments 
|3|-[5|] have so far provided only upper limits on the rate of this decay, recent limits || 
combined with other inputs on neutrino physics might already lead to important information 
on whether neutrinos are Dirac or Majorana particles. These other inputs are indications 
from the analysis of the cosmic microwave background for the presence of a hot dark matter 
component which is presumably neutrinos of mass ~ 1 eV |],|§|| and the indications [|10| ^T~2 



from the analysis of solar and atmospheric neutrinos that neutrinos do oscillate and that 
the mass differences among the three neutrinos are much smaller than this scale of 1 eV. 

We show in this note that if neutrinos are Majorana particles, a combined study of all 
the above pieces of data leads to rather stringent restrictions on two of the mixing angles 
that occur for three flavors of neutrinos. If these results are then confronted with the values 
of these mixing angles allowed by solar, atmospheric and reactor neutrino data, the allowed 
regions are further narrowed and in fact, in a few cases, one is already close to contradiction, 
thus leading to the conclusion that neutrinos are not Majorana particles. 

In the three flavor mixing scheme the neutrino flavor eigenstates v a = z/ e /JjT are related 
to the mass eigenstates = ^1,2,3 by 

Va = ^U ai Vi (1) 

i 

where U a i are the elements of the unitary mixing matrix U. We note that for Majorana 
neutrino [0,[T3] there are three CP violating phases in contrast to the case of Dirac neutrino 
which has only one phase. We use the parametrization 



U 



\ S^S^c' ^ C^jC^S QjS-^j6 ' 3 S^C^pS (jjC^ ^ CipC-(f) J 



(2) 



where 81, 82 and 63 are the three CP- violating phases and c and s stand for sine and cosine 
of the associated mixing angle u, <fi or ip placed as subscript. 

The rate for the neutrinoless double beta decay depends on the following combination 
of the neutrino parameters fl2j: 



E ViU, 



i=l,2,3 



(3) 



where m Vi are the Majorana neutrino masses, U e i the elements of the first row of the mixing 
matrix given in equation (2), rji = (1/i) r\( p = ±1 and r]f p is the CP parity of Majorana 
neutrino z/j. Although neutrinoless double beta decay has not yet been seen experimentally, 
the experimental upper limits on this rate have recently improved to a significant extent. 
In particular one may refer to the results of the Tellurium || and Germanium experiments 
The strongest upper limit so far comes from the Germanium experiment ||5|,|14| and it 

is 



(4) 



m ouf3p < 0.56 eV (99% confidence level) 
< 0.46 eV (90% confidence level) 



These numbers have been obtained using the nuclear matrix elements calculated in [15 



We shall take into account the uncertainties in this calculation (see below). Using the 
conservative number 0.56 eV we have 



E Vi U e 



m. 



i=l,2,3 



< 0.56 eV 



(5) 



Next we consider the fits to the recent data on the anisotropies of the cosmic microwave 
background radiation [16|| and the large scale structure of the universe Best fit |)]|| 

requires a mixture of 10% ordinary baryonic matter, 70% cold and 20% hot dark matter 
with Q m = 1. If the hot component is identified with neutrinos , the model implies [[| 



E 

i=l,2,3 



m 



5eV. 



(6) 



Right hand side of (6) is not expected to be less than 3 eV for Q m = 1 [0] as otherwise 
there is too much small scale power ||18||. On the other hand, solar and atmospheric neutrino 



data suggest that the two mass-squared differences among the three neutrinos are very small 
p0|-p]2[ : m\—m\^ 10~ 5 eV 2 or smaller and m 2 — m 2 ^ 1CT 3 — !CT 2 eV 2 . Hence we take all 



three neutrinos as almost degenerate in mass and using (6), 

m u . « m u w 1.7 eV (7) 

We shall allow m v to vary over a range around 1.7 eV. This will take care of the uncertainties 
of cosmological models as well as those of the calculations of the nuclear matrix elements 
in double beta decay, since only the ratio 0.56/m^ enters into our analysis. Any possible 
improvement in the neutrinoless double beta decay limit can also be incorporated by scaling 
m u appropriately. 

Combining all the inputs, we have the basic inequality 

(m cos 2 u + m sin 2 u e~ i2Sl ) cos 2 + n 3 sin 2 e~ i2S2 < — — (8) 
v ' ' vn v 

where m v is expressed in eV. One can rewrite this inequality in terms of two effective 
phases by combining r\i and Si. However, to make our discussion on CP conserving and CP 
violating cases more transparent we have kept both rji and Si above. We proceed to extract 
the bounds on the mixing angles u and <fi implied by this inequality for various choices of rji 
, Si and typical values of m u favored by the cosmological models. It is to be noted that in 
contrast to the usual oscillation phenomena studied in neutrino physics, CP violation plays 
an important role in neutrinoless double beta decay. 

We shall first consider Si = S2 = in (8). Out of eight possible combinations for different 
values of rji in (8), four combinations are equivalent to the other four, as only the overall 
magnitude in the left hand side of this inequality is constrained. So we shall analyse (8) on 
the basis of four cases : Case I : rji — r\i — 773 = ±1 ; Case II : r\i — —r] 2 = r/ 3 = ±1; Case III 
: —rji = 7/2 = //3 = ±1 and Case IV: T]\—r\i — —7/3 = ±1. Case I is the natural choice for 
rji, if there exists a symmetry linking the three generations. But then the left hand side of 
(8) is unity and so unless m v < 0.56 eV, the inequality cannot be satisfied PD|. Since such 



low values of m„ are not expected in the cosmological models, we conclude that the case 



of equal intrinsic CP parities for the three Majorana neutrinos is not favored. The allowed 
regions in uj and <p for cases II, III and IV are plotted in figures 1, 2 and 3 respectively for 
m u = 1.7 eV. For both II and III, small <fi (<p ~ 0) requires u to be in the region of 45°. 
Whereas small to in II is not allowed at all, small to in III requires large 0. Case IV leads to 
a constraint condition on 0, that is independent of uj and requires to be in the region of 
45°. Figures 4 and 5 show the total allowed regions for all possible combinations of rji (i.e., 
all the cases I, II, III and IV) for m v = 6 eV and 0.64 eV respectively. One may note that 
for movpfi/rriv — > 0, the allowed values of uj and are constrained to lie on the three curves 
: tan 2 = ±cos2u;, 1. Our figure 4 (where mQ V pp/m v = 0.093) gives a small width to these 
curves. 

We next consider CP violating case. Now, the choice rji = r\i = 773 = ±1 , is not ruled 
out. Thus CP violation is capable of changing the conclusions dramatically. Figures 6-8 
show the allowed regions in uj and for a few choices of the parameters 81 and 62- 

We have not considered "maximal" CP violation 5\ = ir/2 and/or 5 2 = tt/2 since these 
cases are subsumed by the choices of relative negative rji, as far as our inequality (8) is 
concerned. Therefore, it is important to keep in mind that our figures (1-5) are relevant 
for such "maximal" CP-violating cases also. In particular, figure 1 and figure 2 can be 
considered to be the "maximally" CP- violating case of Si = 7r/2 and 62 = 0. 

Let us now compare the above bounds on the mixing angles with the results of the 
analysis of data from reactor, solar and atmospheric neutrino experiments. Before we make 
this comparison, we must comment on the role of CP violation in these experiments. We 
note that even for Majorana neutrinos, the oscillation phenomena are controlled by a single 
CP- violating phase 6 pl| . 



The CHOOZ reactor experiment |22] interpreted within a three flavor framework [23 
leads to the constraint 



< 12.5° 

Although CP violation was neglected in reference 
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(9) 



, it is easy to see that (9) is valid 



even if CP is violated. Comparing (9) with figure 3, we see that case IV is ruled out, while 
comparisons with figures 1 and 2 rule out all values of uj except those in the region of 45°. 
These are for m v = 1.7 eV. The results are weakened for m v = 0.64 eV (figure 5) and are 
further strengthened for m v = 6 eV (figure 4). If CP is violated, with 5± = and 62 7^ 
(figure 6), the above conclusion is still valid. But for m u = 1.7 eV, 7r/8 ~ 5i ~ 3n/8 and 
§2 = 0, ( figure 7) as well as n/8 ~ 5\, 82 ~ 3ir/8 (figure 8) there is no allowed region at all 
after using (9). 

It is important to note that as a consequence of (9), the effect of CP violation in all 
the neutrino-oscillation phenomena is much reduced since the CP violating phase factors 
e ±lS always occur in combination with sin0, as sin0 e ±lS . Hence it is legitimate to compare 
with the results of analyses of solar and atmospheric neutrinos even though CP violation is 
usually ignored in those analyses. 



Analysis of solar neutrino data on the basis of the MSW effect leads to the allowed 
uj values of either 2° — 3° or 20° — 40° for <fi < 12.5°. So comparing with the above result, 
we conclude that the result of the present analysis contradicts the small uj MSW solution 
(uj ?z 2° — 3° ), but there exists some overlap with the large uj MSW solution (uj « 20° — 40°). 
Also there is no contradiction with the vaccuum oscillation as a solution of the solar neutrino 
problem since this also requires large uj . Finally there is no contradiction with the results 



of the atmospheric neutrino analysis |24] since this does not involve uj 



So far we have used the information on m u from cosmology to get results on the mix- 
ing angles which were then compared with the results of reactor, solar and atmospheric 
neutrinos. In view of the uncertainties of cosmological models, one can ask what kind of 
information on the quasi-degenerate mass for Majorana neutrinos can be obtained from our 
analysis, if we drop the cosmological input completely. It is clear from our figure 5, that 
for the small uj MSW solution of the solar neutrino problem, we get an upper bound on m v 
of about 0.7 eV and we have checked that this upper bound becomes 4 eV for the large uj 
MSW solution. 

The quantitative results of our analysis are contained in figures 1-8 in the form of re- 



strictions in lo and 0. We may also state two qualitative conclusions that emerge from our 
analysis. 

1. If neutrinos are Majorana fermions and CP is conserved, all the three neutrinos cannot 
have same CP parities. This conclusion (which perhaps is well-known and is included here 
only for the sake of completeness) may be important for model building. 

2. If neutrinos are Majorana fermions, the mixing angle uo cannot be small. Hence, if 
the small uj solution turns out to be the only correct solution of the solar neutrino problem, 
then neutrinos cannot be Majorana fermions. This will have serious consequences for models 
intended to explain small neutrino masses. 

The above results and conclusions are based on the present indications that the exper- 
iments on the neutrinoless double beta decay require wiq^ to be less than a fraction of 
an eV while models with neutrinos as the hot component of dark matter require m v to be 
higher than about 1 eV. The restrictions become more severe and the conclusions become 



stronger if the upper limit on movpp decreases |26| and/or m u increases. 

After the first submission of our manuscript, we came across the papers of Georgi and 
Glashow |28[ and of Branco, Rebelo and Silva-Marcos whose contents have partial 



overlap with our work. Our analysis is more general than both these works in which is 
put zero and further we have considered the CP- violating cases in detail. 
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FIG. 1. The allowed region in (uj,<fi) is shown shaded for m v = 1.7 eV, 
and S± = 62 = 0. 
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FIG. 2. Same as figure 1 but —771 = 772 = % = ±1 
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FIG. 3. Same as figure 1 but 771 = 772 
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FIG. 4. Allowed region for all possible combinations of rji, m v = 6 eV, and <5i = 82 = 0. 
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FIG. 5. Same as figure 4 but for m u = 0.64 eV. 
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FIG. 6. Allowed region for all possible rji, m v = 1.7 eV and 8\ = 0,7r/8 < 82 < 37r/8. 
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FIG. 7. Allowed region for all possible rfi, m v = 1.7 eV and tt/8 < 5\ < 37r/8, 62 = 0. 



10 20 30 40 50 60 70 80 90 



FIG. 8. Allowed region for all possible r/j, m u = 1.7 eV and tt/8 < 61,62 < 3tt/8. 
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